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The electroweak gauge symmetry is allowed to be spontaneously broken by the strongly interact- 
ing vector-like gauge dynamics. When the gauge coupling of a theory runs slowly in a wide range 
of energy scale, the theory is extremely interesting. This may open up the possibility that the 
origin of all masses may be traced back to the gauge theory. We use the Schroedinger functional 
method to determine the scale dependence of the gauge coupling of 10-flavor QCD. Preliminary 
results are reported. 
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1. Introduction 

The main goal of Large Hadron Collider (LHC) is to confirm the Higgs mechanism and to 
find particle contents and the physics low above the electroweak scale. So far many new physics 
models beyond the standard model have been proposed. Among them, Technicolor (TC) [p is 
one of the most attractive candidates \^ as it does not require elementary scalar particles which 
cause, so-called, the fine-tuning problem. This model is basically a QCD-like, strongly interacting 
vector-like gauge theory. Therefore, lattice gauge theory provides the best way to study this class 
of model[^], and the predictions can be as precise as those for QCD, in principle. 

The simple, QCD-like TC model, i.e. an SU(3) gauge theory with two or three flavors of 
techniquarks, has been already ruled out by, for instance, the S-parmeter[Q] and the FCNC con- 
straints. However, it has been argued that, if the gauge coupling runs very slowly ("walks") in 
a wide range of energy scale before spontaneous chiral symmetry breaking occurs, at least, the 
FCNC problem may disappear [Q]. Such TC models are called walking technicolor (WTC) and 
several explicit candidates are discussed in semi-quantitative manner in [g]. Since the dynamics in 
WTC might be completely different from that in QCD and hence the use of the naive scaling in 
N c or Nf to estimate various quantities may not work, the 5-parameter must be evaluated from the 
first principles Although really important quantity is the anomalous dimension of xj/iff opera- 
tor, to find theories showing the walking behavior is a good starting point. Recently many groups 
started quantitative studies using lattice technique to answer the question what gauge theory shows 
walking behavior. In [H], the running couplings of 8- and 12-flavor QCD are studied on the lattice 
using the Schroedinger functional (SF) scheme[Q]. Their conclusion is that while 8-flavor QCD 
does not show walking behavior 12-flavor QCD reaches an infrared fixed point (IRFP) at g 2 R ~ 5. 
In spite of the scheme-dependence of running and its value of IRFP, the speculation inferred from 



Schwinger-Dyson equation [|T0[] suggests that gj R ~ 5 is too small to trigger spontaneous chiral 



symmetry breaking. Although 12-flavor QCD is still an attractive candidate and is open to de- 



bate IJllJl, we explore other Nf. In the following, we report the preliminary results on the study 
of the running coupling in 10-flavor QCD. Since the lattice conference, statistics is increased by 
much. The following analysis is based on the increased statistics. 

2. Perturbative analysis 

Before going into the simulation details, let us discuss some results from perturbative analysis. 
In this work, we adopt the /3 function defined by 

f3(g 2 (L)) =L^^ = b lg \L)+b 2 g 6 (L) + b 3 g\L)+b 4 g l0 (L) + .-., (2.1) 
where L denotes a length scale. The first two coefficients are scheme independent, and given by 



(4n) 2 



11 
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(2.2) 



The higher order coefficients are scheme-dependent and are known only in the limited schemes. In 
this section, we analyze the perturbative running in the following four different schemes/approximations: 
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Figure 1: L-dependence of g 2 (L) in different approximations or schemes. Nf dependence is also shown. 



i) two-loop (universal), ii) three-loop in the MS scheme, iii) four-loop in the MS scheme, iv) three- 
loop in the Schroedinger functional scheme. The perturbative coefficients relevant to the following 
analysis are 



-MS = ^ 

3 " (4k) 6 



2857 5033 325 2 



18 3 54 



(2.3) 



■ MS - 2 [ 29243.0 -6946.30N/ +405.089^+1.49931^], (2.4) 



h SF _ ,MS i b 2 c 2 b l ( c 3 ~ C f) n 

where the coefficients c\ and c\ depend on the spatial boundary condition of the SF used in calcu- 
lations, i.e 6. Those for 6 = n/5 and cf for 6 = are known as 

C2 /5 = 1.25563 + 0.039863 xN f , (2.6) 
c* /5 = ( C 2 /5 ) 2 + 1.197(10) + 0.140(6) xN/- 0.0330(2) xNj, (2.7) 
c° 2 = 1.25563 + 0.022504 x N f , (2.8) 

but C3 is not. Therefore, the case iv) is studied with 6=n/5. It should be noted that in our numerical 
simulation 6=0 and thus, rigorously speaking, the example iv) is not applied to our numerical result. 
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N f 


4 
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10 


12 


14 


16 


2-loop universal 








27.74 


9.47 


3.49 


0.52 


3-loop SF 


43.36 


23.75 


15.52 


9.45 


5.18 


2.43 


0.47 


3-loop MS 




159.92 


18.40 


9.60 


5.46 


2.70 


0.50 


4-loop MS 






19.47 


10.24 


5.91 


2.81 


0.50 



Table 1: The IRFP from perturbative analysis. 



Fig. |] shows L-dependence of g 2 (L) in the four different schemes/approximations, i)-iv), 
with the initial condition g 2 (Lo) = 1.6. The values of perturbative infrared fixed point (IRFP) 
are summarized in Tab. [T[ As seen from Tab. once one goes beyond the two-loop approxi- 
mation, the fixed point value of 12-flavor QCD is stable at around g 2 R ~ 5 against the change 
of schemes/approximations. It is interesting that this IRFP is completely consistent with that 
obtained in the non-pertrurbative calculation by Now looking at the perturbative IRFP at 
Nf = 10, the similar stability is seen at g 2 ~ 10. Furthermore, according to an analysis based on 
the Shwinger-Dyson equation, there is an argument that chiral symmetry breaking occurs at around 
g 2 ~ An 2 1 (3C2(R)) = 7Z 2 JTO]]. In summary, the perturbative analysis suggests that 10-flavor QCD 



is the most attractive candidate for WTC. 



3. Simulation details 

We employ the Schroedinger functional method [§] to determine the running coupling con- 
stant. Unimproved Wilson fermion action and the plaquette gauge action without any boundary 
counter terms are adopted to describe regularized dynamics of techniquarks and technigluons, re- 
spectively. The parameter of the spatial boundary condition for fermions, 6, is set to 0. To de- 
termine the scale dependence, calculations on several different lattice sizes are required. In this 
analysis, we report the results obtained with (L/a) 4 = 4 4 , 6 4 , 8 4 and 12 4 lattices. The calculation 
on 16 4 lattice is in progress. The bare gauge coupling j8 = 6/g^ is explored in the range of 4.4- 
24.0. At around /3 ~ 4.4, we encounter a bulk phase transition, where the plaquette value suddenly 
jumps to a smaller value. Whenever this happens, we discard the configuration. The numerical 
simulation is carried out on several architectures including GPGPU and PC cluster. The standard 
HMC algorithm is used with some improvements in the solver part. Since the Wilson type fermion 
explicitly violates chiral symmetry, the critical value of k has to be tuned to the massless limit. We 
performed this tuning for every pair of (j3, L/a). So far, we have accumulated 5,000 to 200,000 
trajectories depending on (j8, L/a). 

The procedure to determine the running is standard and consists of the following steps. 

1. Calculate the Schroedinger functional scheme coupling at various g^'s on several L/a's, 
which is denoted by g|/ a (,?o)- 

2. Fit g 2 L i a {go) at each L/a as a function of g^ to obtain an interpolating formula g 2 ^(gl)- 

3. Set the input value uo, which implicitly sets an initial length scale Lo- 

4. Choose one L/a and determine g 2 , * satisfying g^Q?o'*) = M o- 

5. Choose one L' /a = sL/a (typically s =2 or 3/2) and read the value of g 2 ^j a {g^ '*)> which is 
stored in L(sLo/a, Lq/o). 
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Figure 2: g^-dependence of gQ/g 2 (L) for L/a = 4, 6, 8, 12. The right panel is just an enlargement of the 
left. 

6. Make a perturbative correction via r(5Lo /a, Lo/a) = g^ a /[l + 8(a/sLo, a/Lo) x«] with the 
known coefficient 8. 

7. Repeat 4-6 with different a, i.e. for L/a' and L' /a' = sL/a! . 

8. Take the continuum limit by using L(sLo/a, Lq/o) and L(sLo/a' , Lo/a') to obtain <j(sL ). 

9. Set u\ = a(sL ) 

10. Repeat 4-9 several times with the initial value u\, U2, ■■■■ Then we obtain a series of 
o(sLq), o(s 2 Lq), a(s 3 Lo), This scale dependence describes the running of the cou- 
pling. 

4. Results 



L/a 


N 


£ 2 /dof 


a L/a,l 


«L/a,2 


a L/a,3 


a L/a,4 


4 


2 


45.6(2.2) 


0.4903( 0.0006) 


-0.4677( 0.0025) 






4 


3 


11.0(1.1) 


0.5003( 0.0034) 


-0.3909( 0.0354) 


-0.1038( 0.0435) 




4 


4 


2.2(0.5) 


0.5080( 0.0008) 


-0.4746( 0.0094) 


0.1168( 0.0232) 


-0.1486( 0.0148) 


6 


2 


19.5(1.6) 


0.4607( 0.0016) 


-0.3826( 0.0038) 






6 


3 


4.3(0.8) 


0.4888( 0.0015) 


-0.2662( 0.0085) 


-0.1856( 0.0115) 




6 


4 


2.1(0.6) 


0.5028( 0.0026) 


-0.4675( 0.0422) 


0.2486( 0.0894) 


-0.2593( 0.0516) 


8 


2 


12.4(1.5) 


0.4822( 0.0018) 


-0.4004( 0.0069) 






8 


3 


1.9(0.6) 


0.5088( 0.0021) 


-0.2135( 0.0192) 


-0.2668( 0.0239) 




8 


4 


0.7(0.4) 


0.5153( 0.0020) 


-0.3610( 0.0515) 


0.0709( 0.1108) 


-0.1972( 0.0605) 


12 


2 


9.7(1.9) 


0.482 1( 0.0047) 


-0.3870( 0.0162) 






12 


3 


2.7(1.0) 


0.5164( 0.0040) 


-0.1716( 0.0316) 


-0.3104( 0.0363) 




12 


4 


1.5(0.8) 


0.5269( 0.0033) 


-0.3959( 0.0740) 


0.2611(0.1728) 


-0.3495(0.1003) 



Table 2: The coefficients determined in the fit. 



Fig. U essentially shows the Schroedinger functional coupling obtained on the four lattices 
and their fit curves. Notice that, instead of g£/ a , we deal with l/g^/ a just for convenience. As an 
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Figure 3: Running of Schroedinger functional coupling. 



interpolating function, we chose 

1 -aL/ a> i8o 



J2n ' 



S L /a(S 2 ) 1 +PWa X^ + I„=2 a L/a,n * 

where Pi L / a is the L/a-dependent, one-loop coefficient, the values of which for various L/a are 



(4.1) 



Pi, L/a 



0.4225003137 for L/a = 4 

0.4477107831 for L/a = 6 

0.4624813408 for L/a = 8 

0.4756888260 for L/a = 12 



We optimize the degree of polynomial in the denominator of (11), N, by monitoring # 2 /dof. The 
results of the fits are tabulated in Tabs. ||. From the # 2 /dof values, N=4 is chosen in all the cases. 

Setting the initial value of the running uq to 3.0, we calculated the running through the proce- 
dure explained above. Since we did not implement any 0(a) improvements, large scaling violation 
is expected. This can be improved only by performing simulations at larger lattices. We tried to 
evaluate the size of the scaling violation, but it fails because of the large statistical error. Thus, 
in this report, instead of the continuum limit, the step scaling function constructed from the data 
on L/a=% and 12 lattices is taken as the current estimate. A subtlety occurs in the perturbative 
correction, c.f. the step 6 in the procedure. While in the weak coupling region o(L) < 2 it works, 
in the other region the correction turns out to enhance the scaling violation. This happens when 
the direction of the scaling violation at the one-loop level is opposite to that of the real violation. 
Thus, we calculate the running coupling with and without the perturbative correction and interpret 
the difference as the potential systematic uncertainty associated with the scaling violation. Fig. || 
shows the result thus obtained. The perturbative predictions with 6 = n/5 is also shown for a 
reference. Although the statistical error is still large, an interesting observation is that the data 
without the perturbative correction almost goes along the 3-loop curves and at g| p ~ 10 turns to 
rapid increase. 
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5. Summary and Outlook 

The running coupling constant of 10-fiavor QCD is studied. The perturbative analysis suggests 
that this theory is extremely interesting. The preliminary result obtained without continuum limit 
shows an interesting behavior. In order to draw definite conclusions, we clearly need larger lattices 
and more statistics. The calculation is in progress. 

Numerical simulations are performed on Hitachi SRI 1000 at High Energy Accelerator Re- 
search Organization (KEK) under a support of its Large Scale Simulation Program (No. 09-05), 
on GCOE (Quest for Fundamental Principles in the Universe) cluster system at Nagoya University 
and on the INSAM (Institute for Numerical Simulations and Applied Mathematics) GPU cluster at 
Hiroshima University. This work is supported in part by the Grant-in-Aid for Scientific Research 
of the Japanese Ministry of Education, Culture, Sports, Science and Technology (Nos. 20105001, 
20105002, 20105005, 21684013, 20540261 and 20740139), and by US DOE grant #DE-FG02- 
92ER40699. 
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